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ABSTRACT 

We examine correlations between the masses, sizes, and star formation histories 
for a large sample of low-redshift early-type galaxies, using a simple suite of dynamical 
and stellar populations models. We confirm an anti-correlation between size and stellar 
age, and survey for trends with the central content of dark matter (DM) . An average 
relation between central DM density and galaxy size of (pdm) oc R e e~ 2 provides the 
first clear indication of cuspy DM haloes in these galaxies — akin to standard ACDM 
haloes that have undergone adiabatic contraction. The DM density scales with galaxy 
mass as expected, deviating from suggestions of a universal halo profile for dwarf and 
late- type galaxies. 

We introduce a new fundamental constraint on galaxy formation by finding that 
the central DM fraction decreases with stellar age. This result is only partially ex- 
plained by the size-age dependencies, and the residual trend is in the opposite di- 
rection to basic DM halo expectations. Therefore we suggest that there may be a 
connection between age and halo contraction, and that galaxies forming earlier had 
stronger baryonic feedback which expanded their haloes, or else lumpier baryonic ac- 
cretion that avoided halo contraction. An alternative explanation is a lighter initial 
mass function for older stellar populations. 

Key words: dark matter - galaxies : evolution - galaxies : general - galaxies : 
elliptical and lenticular. 



1 INTRODUCTION 

The formational history of early-type galaxies (ETGs: el- 
lipticals and lenticulars) remains an outstanding question. 
While these dynamically hot systems may be basically un- 
derstood as end-products of galaxy mergers, the details 
of these mergers and their cosmological context are un- 
clear. High-redshift (z) observations have made initially sur- 
prising discoveries that many ETGs were already present 
at early times with mature stellar populations, and that 
these galaxies we re much more compact than those in the 
present day, (e.g. iGlazebrook et al.ll2004l : iDaddi et al1l2005l : 
iTruiillo et al.ll2006l ). 
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The evolution in ETG sizes is st ill controversial in 
both observation and interpretation (e.g. Ivan Dokkum et al.l 
120091 . l20ld ; iMancini etail l20ld : IValentinuzzi et alj I2OI0F T 

However, the most likely scenario is for a combination of 
effects where individual galaxies grow in size by accretion 
of smaller, gas-poor galaxies (an "inside-out" picture of 
galaxy formation), and where younger ETGs are formed 
with larger sizes because of their decreased cold gas content 
and the lower background densities of dark matter (DM; e.g 



Saracco et al.ll2009l:lBezanson et "al]|2009l: Ivan der Wei et al 
20091; iNaab et al J 120091: iHopkins et alj|20ld: IShankar et a! 



2010) 



The role of DM is considered fundamental to the for- 
mation of galaxies, and DM halo properties have been 
extensively studied in cases such as gas-rich spirals and 
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near by dwarfs which have suitable observa t ional tracers 
(e.g.lMcGaugh et alj|2007l ; IWalker et alll2009l ; iKalirai et all 
120101 ). Studying DM in the general population of ETGs is 
in many ways more difficult, with ongoing surveys of the 
available large-radius halo trace rs attempting to remedy 
our ignorance in thi s area (e.g. iRqmanowskv et alJ 120091: 
Coccato et al] 120091 ; IProctor et al.l |2009|; iThomas et ail 



20091 . hereafter T+09; IWoodlev et al.ll201of ). 

DM can alternatively be studied with less precision 
but in more extensive ETG samples by considering the 
well-studied central regions — inside the "effective radius" 
(Reg) enclosing half the stellar light, where DM is generally 
thought to be a minor yet potentially detectable contrib- 
utor to the mass. Here one of the classic approaches is to 
analyze the "fundamental plane" (FP) relating ETG sizes, 
luminosities, and central velocity dispersions (00). The FP 
shows a "tilt" or systematic deviation from simple expec- 
tations based on galaxies with constant dynamical mass-to- 
light ratios M/L, probably implying systematic differences 
in the stellar populations or in DM content. 

After many years of debate, there i s still not a consen- 
sus on what is driv i ng the FP tilt (e.g. iTruiillo et al.1 12004 
Ferreras et al.ll2005l;ICappellari et alj|2006l. C+06 hereafter; 



Bovlan-Kolchm et alJIgOOAlDekel fc Coxj|2q06l:lBolton et all 
2001 l2008l; lJun fe Iml 120081; iTortora et al.l 120091. hereafter 
pape r I; lAllanson et al.ll2009l; iGraves! 120091; (Grillo fc Gobatl 



bold ; iHumphrev fc Buotd l20ld ; iLa Barbera et al.1 l2010al ) 



Some work suggests stellar populations variations or non- 
homologies in the luminosity profiles as the tilt drivers. How- 
ever, it is a fairly generic expectation from the standard 
cosmological framework for galaxy formation that the cen- 
tral DM content of ETGs will systematically increase with 
luminosity — a point we discussed in paper I and develop fur- 
ther in this paper. If the FP tilt is not caused in large part by 
DM, there could be problems implied for galaxy formation 
theory. Here one could pursue two different philosophies: to 
empirically and robustly determine the reasons for the FP 
tilt without recourse to theory (e.g. making no assumptions 
about the underlying DM profiles); or to adopt the theoret- 
ical framework as broadly correct and consider the detailed 
implications for ETG composition and formation. 

Following the second approach, it is now time to be- 
gin moving on from phenomenological questions about the 
FP tilt, and to establish more direct connections between 
DM in ETGs and their formational histories. In particu- 
lar, the central DM content could prove crucial to solving 
the size-evolution puzzles mentioned above, and more fun- 
damentally to understanding the assembly of ETGs. To this 
end, we now extend the analysis of paper I, which combined 
models of stellar dynamics and population synthesis to infer 
total and stellar masses in a large sample in nearby galaxies, 
and thereby to analyze the FP tilt. There have been previ- 
ous suggestions that age and star form ation timescales are 
impo r tant fourth parameters in the FP ( Terlevich fc Forbes! 
120021 ; iGargiulo et al.ll2009l ; IGraves! 120091 ) ■ We will now ex- 
plore this possibility systematically, using the results from 
paper I to consider additional correlations involving DM 
content and star formation histories (SFHs). 

It should be noted at the outset that the data and the 
analysis techniques that we use for deriving mass and SFH 
parameters may not be the most state-of-the-art. However, 
our aim is to pioneer a framework for interpreting any data 



of this kind in a broad cosmological context (where in par- 
ticular, emerging high- z data-sets provide only crude obser- 
vational constraints on ETGs), and we are so far able to 
tentatively identify some basic and intriguing trends. 

The paper is organized as follows. In Section [2] we 
present our basic observational results. We go on to ana- 
lyze some implications of the trends of DM with mass in 
Section [3] and with age in Section [4] In Section [5] we sum- 
marize our conclusions. Two Appendices include analyses of 
statistical and systematic uncertainties. 



2 OBSERVATIONAL RESULTS 

Here we present our basic observational results. Section \'2.1 1 
provides an overview of our galaxy sample and mass infer- 
ences, and summarizes the trends versus mass. Section [2.21 
presents results related to galaxy age. 



2.1 Summary of sample and initial results 

Our starting point is a c ollection of 335 local ETGs from 
IPrugniel fc Simie 3 <| 19961 ) that we recently re-analyzed in 
paper I. The sample was selected to have measurements of 
(To and at least two colours, with a subsample of ~ 220 
systems that also include the maximum rotation velocity 
(Knax). There are 218 elliptical galaxies and 117 lenticu- 
lar/SO systems: as shown in paper I, the two subsamples 
show similar stellar and DM properties and we will consider 
them jointly in this work. 

We next summarize the main steps of the paper I anal- 
ysis, starting with the stellar population models. We used 
a set of simple stellar po pulation (SSP) syn t hetic spectra 
from the prescription of iBruzual fc Chariot! (|2003l . BC03 
hereafter ) to fi t th e observed galaxy colours. W e assumed a 
ISalpeteij j 19551 ) or IChabriej (|200ll .' l2002l . 120031 ) initial mass 
function (IMF) alternatively, with initial masses m in the 
range 0.1 — lOO Mp). For the c urrent paper, we instead use 
an intermediate iKroupal l|200ll ) IMF as our default mode{]]. 

These single-burst models were convolved with an expo- 
nential SFR oc e~ t//T to generate more general SFHs, where 
t is a characteristic time scale. The age, metallicity (Z), 
and r were free parameters of the model while stellar M/L, 
T t , was inferred from the best-fit solution for each galaxy. 
The reliability of the modelling technique and the intrinsic 
parameter scatter, as well as the presence of spurious cor- 
relations among the stellar parameters induced by the stel- 
lar modelling procedure, have been checked through Monte 
Carlo simulations (see paper I for details, as well as Ap- 
pendix [B3] of this paper). In addition, a recent novel, com- 
pletely independent technique for estimating T* using glob- 
ular cluste r systems has yiel ded results in perfect agreement 
with ours (|Forte et al.ll2009l . Fig. 6). 

We derived the dynamical M/L (Td yn ) within _R e ff 
by means of Jeans analysis assuming spherical symmetry, 
isotropy of the velocity dispersion tensor, and introducing a 



For old stellar populations ( J>1 Gyr), changing the IMF gen- 
erally impacts the colours at less than the ~ 0.01 mag level. The 
IMF change can then be fairly represented as a simple overall 
stellar M/L renormalization, which is reduced by a factor of 1.6 
in changing from Salpeter to Kroupa; see e.g. Fig 4 of BC03. 



Early-type galaxy formation and dark matter 3 



25 



20 



15 



o 



10 



. " : \* *•'• 
• t. . •. • •• 



_L 



1 1 



4 

T./T 



Figure 1. Stellar versus dynamical M/L within R c g for the early- 
type galaxy sample as derived in paper I, using a Kroupa IMF. 
The dotted line is the one-to-one relation, and typical statistical 
uncertainties are illustrated by the error bars at left. For most of 
the galaxies, the presence of a significant mass excess is evident, 
implying a non-zero DM fraction in the central i? e ff- 



rotation velocity correction in the spherically averaged RMS 
velocity drms = Vv 2 + a 2 , which is a measure of the total 
kinetic energy of the system within R c g. In the following 
we will use Td v n obtained w ith the multi-component model 
in paper I (i.e. ISersid [l968l profile for the stars and singu- 
lar isothermal sphere for the total mass); we have verified 
that the particular mass model choice does not affect the 
conclusions of this paper. 

Central DM fractions are inferred by stellar and dynam- 
ical M/L estimates, being by definition 



/dm 



Mi, 



M~o 



(1) 



Mtot Mtot Tdyn 

where and Td yll are obtained in paper Q the DM mass 
includes any mass in diffuse gas. Note that as seen in Fig. (TJ 
a fraction of the galaxies (under any of our IMF choices) ap- 
pear to have /dm < 0, which would be an unphysical situa- 
tion. This is not necessarily a worry, since errors in the M/L 
estimates can scatter some data to /dm < 0. We have in- 
vestigated this issue quantitatively in AppendixfX] finding a 
strong but not definitive suggestion that some ETGs are not 
compatible with having a Salpeter IM F. Simi lar points were 
made in other studies (|Renzinill2005l ; C+06; iFerreras et al.l 
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Figure 2. Relation between central DM fraction /dm an d stel- 
lar mass Mi, for ETGs, for an assumed Kroupa IMF. Individual 
galaxy data points are shown, colourised according to their ages 
as illustrated by the colour-bar at lower left; the local age varia- 
tions in the plots are smoothed in order to see the general trends. 
Except for the lowest masses (which are dominated by SOs), /dm 
increases on average with M* (illustrated by the black points with 
error bars showing the binned averages with scatter). The curves 
show various predictions of ACDM toy models (see Section 13.11 1. 
The bottom set of 9 curves shows model predictions without adi- 
abatic contraction (AC), for three age bins (3, 7, 11 Gyr, again 
colourised by age), and three egF values (0.7, 0.3, 0.1: solid, long- 
dashed, short-dashed respectively). The top set of curves show 
predictions with AC, including a smaller variation of model pa- 
rameters for the sake of clarity. Long-dashed curves show the three 
age bins again, with esF=0.3 and the G+04 AC recipe. Solid and 
short-dashed curves show the egp=0.7 and 0.1 cases for G+04 
AC and 7 Gyr. The heavy long-dashed curves show egp=0.1 for 
B+86 and 7 and 3 Gyr. 



l2008l : ICappellari et al.ll2009l : iBarnabe et al.l feoiO'l using con- 
ceptually similar techniques. Treu et all (|2010l . hereafter 
T+10) on the other hand claimed from an analysis of gravi- 
tational lenses that a Salpeter IMF is preferred, with a possi- 
ble systematic IMF variation. We will return to these issues 
in later Sections, and for now note that an assumed universal 
Kroupa IMF is a reasonable starting point. 

We found in paper I that the /dm value correlates with 
luminosity and stellar mass (Fig. [2}, a trend that could be 
the phenomenological cause of the FP til10. We also exam- 
ined the formational causes of the tilt in paper I (Section 
6), using toy models of galaxies in a cosmological context 
to verify that the expected scaling relations of DM haloes 



2 While T(j yn is a three-dimensional quantity, T* is projected 
and susceptible to additional variations from large-radius material 
along the line of sight. However, we have estimated using some 
simple models that T* would typically be affected at only the 
~ 1% level. 



3 In principle, the tilt could instead be caused by a systematic 
variation of IMF with galaxy luminosity (see paper I). We will 
reconsider IMF variations later in this paper. The effects of "non- 
homology" of the stellar profiles were implicitly corrected for in 
our modelling. 
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Figure 3. Galaxy stellar age versus stellar mass. The points are 
colour-coded according to the DM fraction (after smoothing; see 
inset colour bar) . The overall average trend in mass bins is shown 
by points with error bars representing the 1 a scatter. 

might naturally explain the /dm trends. We will develop 
this theme further in the following Sections, and also begin 
studying the implications of the SFHs. 

One basic SFH result is illustrated by Fig. [3] the stel- 
lar ages and masses of the galaxies are correlated, providing 
another example of the well-known "archaeological downsiz- 
ing" phenomenon wherein more massive systems form earli- 
est. We plan to consider the physical reasons for this down- 
sizing in a subsequent paper, comparing cosmologically- 
based models for SFHs and DM halo assemblies. Here we 
will simply treat the SFHs as a given which we will attempt 
to correlate with DM properties, envisioning these as funda- 
mental ingredients for constraining future theoretical mod- 
els. 

2.2 New trends with age 

Here we move beyond the FP analysis of paper I and ex- 
amine in more detail the DM trends, including the connec- 
tions with SFHs. We plotted in Fig. [2] the empirical values 
for /dm versus stellar mass; the two quantities are corre- 
lated albeit with substantial scatter and suggestions of non- 
monotonicity. We have also colour-coded the data points in 
Fig. [2] by stellar agj3, an d can see immediately that much of 
the scatter in /dm might be explained by systematic correla- 
tions with SFHs. In particular, at a fixed mass the galaxies 
with the youngest stars are found to have on average the 
highest DM fractions. This is a central result of our paper 
which we consider in more detail below. Trends involving r 
will be studied in a subsequent paper. 

We show two other projections of the three-dimensional 

4 Hereafter, where applied, the smoothed colour coding is ob- 
tained by linearly interpolating the average values of the third 
(colour-coded) quantity over a regular grid of the other two plot- 
ted quantities. This procedure allows us to readily see qualitative 
trends in the data, with the caveat that quantifying these trends 
can be unreliable after smoothing. 
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Figure 4. Central DM fraction of ETGs versus stellar age. The 
results for a Kroupa IMF are indicated by the left-hand axis, and 
for Salpeter by the right hand. The data points are colour-coded 
by the locally-averaged stellar mass, with colour coding as in the 
inset. A typical uncertainty ellipse for an individual galaxy (at 
the centre of the range) is shown at upper right. 

space of age-mass- /dm in Figs. [3] and 2] The combination of 
/DM-mass and age-mass correlations (Figs. [2] and [3} suggest 
that there could be an overall /oM-age correlation. However, 
Fig-l3]shows that there is a net anta-correlation between /dm 
and age, which implies that /dm couples more strongly to 
age than to mass (we will revisit this issue in Section H. 11 1. 

Before analyzing these trends in more detail, we con- 
sider that there is an important fourth parameter involved, 
the effective radius R c s. This is because /dm is evaluated 
at iicff, a variable parameter that has its own systematic 
dependencies. As we found in paper I and will revisit be- 
low, the well-known ETG size-mass relation (R a a-M t ) can 
explain much of the /dm-M* trends as an "aperture effect" , 
without invoking any intrinsic variation in the DM proper- 
ties. This is because the stellar _R ff varies more quickly with 
mass than do the inner DM halo properties (in ACDM), 
meaning that for a more massive galaxy, the larger R e g en- 
closes a higher fraction of the total DM. 

The next question is whether there is a link between 
R e f! and age that could be driving the /DM-age trends. As 
mentioned in Section [TJ the size-mass relation is thought 
to evolve strongly with time, such that ETGs at a fixed 
mass are more compact at earlier times. This trend can be 
considered generically to include contributions both from 
size growth of existing galaxies, and from the birth of new 
galaxies with systematically larger sizes. As a consequence 
of the latter effect, we would expect the population of ETGs 
at a fixed time (e.g. z = 0) to show a correlation between 
size and age at fixed mass, such that the older galaxies are 
smaller. 

This qualitative prediction has recently been confirmed 
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Figure 5. Size versus age relation for ETGs at z = 0, with galax- 
ies colour-coded by their stellar mass. Solid lines show the average 
trends in these mass bins and are labelled with the log 10 (M^/Mq) 
values (Kroupa IMF). Open and filled symbols show objects that 
have lower and higher /dm (respectively) than the average for 
other galaxies of similar mass and age. 

by low-z observations (SB09; Ivan der Wei et all 12009). and 
we now investigate it for our own ETG sample. We plot R c a 
against age in Fig. [5] where galaxies are colour-coded ac- 
cording to four fiducial stellar mass bins that we will use 
throughout the paper: log 10 (M*/M Q ) ~ 10.3, 10.8, 11.2, 
11.6 for a Kroupa IMF). The overall galaxy distribution does 
not show any significant dependence of i? e ff with age, but 
for any given mass there is a clear anti-correlation (particu- 
larly for the most massive galaxies), i.e. older galaxies have 
on average smaller effective radii. We quantify the size-age 
dependency by fitting a log-linear relation in each of the 
mass bins. The slopes range from —0.005 to —0.030, where 
the normalization increases with mass in accordance with 
the typical R e s-M* relation (see e.g. paper I). These trends 
appear to be independent of the galaxy sub-type (ellipti- 
cals versus SOs) — which is a crucial point because one might 
otherwise speculate that the younger objects are SOs with 
systematically different properties than ellipticals. 

These size-age trends might be driven by higher den- 
sities of DM and gas at earlier epochs of galaxy formation 
(see references in Section [1]). However, as with the age-mass 
correlation, we will for now set aside a physical interpreta- 
tion of the size-age correlation and simply consider it as a 
nuisance factor that must be treated appropriately in our 
later toy models of DM content. These size-age results are 
not meant to be universally applicable to ETGs but are a 
characterization of this particular data-set as needed for our 
self-consistent galaxy toy models. 

The size-age trends clearly go in the right direction 
to explain the /DM-age trends: older galaxies could exhibit 
lower /dm simply because their R c s values are smaller, so 



smaller volumes of their DM haloes are probed. In princi- 
ple, one might be able to see directly from the data whether 
/dm is more directly coupled to age or to ReS, and to this 
end we have highlighted the high- and low- /dm galaxies in 
Fig. [5] However, it is not immediately obvious whether R e g 
or age is the stronger factor, probably because of the scat- 
ter in the data. Instead, in the next section we will adopt 
a model-dependent approach, attempting to make sense of 
the observed trends in the context of standard DM+galaxy 
models. 

Before continuing, in Appendix [B] we carry out a num- 
ber of cross-checks on the robustness of our basic /DM- 
age result. First, comparing independent constraints from 
the literature we find some support for an /DM-age anti- 
correlation. However, there may be disagreement in the 
residual /oM-age trend after accounting for the -R e ff-age 
relations (which we examine further in the following Sec- 
tions). Such differences illustrate the need for further in- 
depth study of this topic. We also check in Appendix [B] 
the effects of our stellar populations modelling assumptions, 
finding that our main conclusions are fairly insensitive to 
these. 

Besides the systematic uncertainties, there is the issue 
of the statistical errors being correlated, in the sense that an 
error in age correlates with an error in T* and produces an 
artificial /DM-age anti-correlation (see error ellipse in Fig. [4}. 
To gauge this effect, we have simulated mock data-sets with 
random errors applied to the ages, propagated these to errors 
on T*, and then regenerated the /DM-age plot. The artificial 
trend generated by the errors turns out to change the slope 
by only ~ —0.005 Gyr -1 (compared to the observed slope 
of ~ —0.04 Gyr -1 ). We will continue in the rest of this 
paper to assume that our empirical results are correct, and 
go on to examine the implications for galaxy structure and 
formation. 



3 DARK MATTER IMPLICATIONS 

Given the preceding inferred DM trends, we begin here 
with interpretations of the mass dependencies, constructing 
ACDM-based toy models in Section 13.11 We consider DM 
density profiles in Section 13.21 and comparisons to related 
literature results in Section [3. 31 

3.1 ACDM halo models and trends with mass 

We now construct a series of toy models of galaxies includ- 
in g DM haloes base d on A CDM theory, as previously done 
in iNapolitano et alj (|2005t ) and in paper I. With this ap- 
proach we are not modelling individual galaxies but gener- 
ating typical representatives in the space of mass and age. 
Each model is spherically symmetric and is comprised of a 
stellar spheroid and a DM halo. The spheroid has a Sersic 
density distribution with an n-index as described in paper 
I, and an i? c ff-M*-age dependence taken from the log-linear 
fits in Fig. 

The DM halo has an NFW density profile 
|Navarro et alj Il997l ) following a typical concentration- 
virial mass (cvir-M v i r ) relatiorQ The DM halo is optionally 

5 As in our previous papers, we have adopted a theoretical rela- 
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"adiabatically contracted", which is an approximate way 
to model the expected drag of dissipatively infalling 
baryons on the surrounding DM, producing a halo with 
a higher central DM density than in collisionless iV-body 
simulations. 

The t wo AC recipes we use a re the most common 
ones, from iBlumenthal et al.l ()l986l . hereafter B+86) and 
iGnedin et all l|2004 hereafter G+04), with the latter pro- 
viding a weaker (and probably more realistic) effect. Other 
recipes are also availab le (e.g. ISellwood fc McGauehl 120051 : 
iGustafsson et all 120061 ). but a more important caveat is 
that the entire slow, smooth-infall AC scenario may not 
be correct for forming cosmological structures — particularly 
the ETGs which may have experienced relatively violent, 
clumpy assembly histories. Current high-resolution simula- 
tions do suggest that the standard AC models are not accu- 
rate in detail, and that the baryonic eff ects on DM haloes 
might be highly var i able and stochastic ( Abadi ct al. 2009; 
Pedrosa et all l2009l . l20ld : iTissera et all l20ld : iDuffv et all 



300 



2010l ). For simplicity, we will consider our AC and no- AC 



models as representing two extreme models of galaxy for- 
mation, where the baryons have alternatively strong or we ak 
net impact on the DM halo fcf. lLackner fc Ostrikcr 2010). 

In order to match the galaxies to their haloes, we must 
relate their stellar masses to their total (or virial) halo 
masses. We can parameterize this connection by assuming 
the cosmological bar yonic fraction, /b ar = f2bar/^m = 0.17 
ijSpereel et all l2007h . is conserved within each halo. The 
stellar mass is then given by an efficiency of star forma- 
tion esF, so that = 0.17esFM v i r . We use a series of 
plausible values e S F = (0.1, 0.3, 7) (e . g. IZheng et al.ll2007l: 
Brown et all 120081; iMore et ail 120091; IConrov fc Wechsler 



20091 ; iMoster et ail Eoidl : iGuo et alJ l20ld : lLagattuta et a! 
20ld : ISchulz et alj|20ld . hereafter S+10). 



We also adopt the Kroupa IMF as our fiducial choice 
but consider the impact of Salpeter or Chabrier IMFs. Our 
model has effectively three free "parameters" : esf, adiabatic 
contraction (AC) recipe and IMF choice. Fig. [6] shows an 
example of a series of galaxy toy models for fixed stellar 
halo masses, and varying AC assumptions. It can be seen 
that the DM content inside 1 R e s varies most by changing 
from a no-AC model to an AC model. The choice between 
the B+86 and G+04 recipes makes only a small difference 
on these scaled- This model may seem simple, but as far as 
we know it has never been applied to the FP before, except 
for the related study of S+10 as we discuss later. 

We now begin by deriving the predicted DM fraction 
within Rtfi as a function of mass, in bins of constant age (i.e. 
fixing the R c g values for a given mass), and with a number 
of different assumptions on esF and AC recipe. The results 
are shown in Fig. [2] it can be seen first of all that without 
AC, the models generally underpredict /dm- The AC models 
on the other hand match the data much better on average, 



tion based on a DM power spectrum normalization of ag = 0.9 
llBullock et alj|200ll) . We have also derived models based on a 
more recent estimate of erg = 0.8 (Maccio et al . 2008). which pre- 
dict lower central DM content, but only at the level of changing 
/dm by ~ -0.05. 

6 The AC recipe differences become more important at smaller 
radii, where the detailed treatment of the stellar mass profile also 
becomes more important, e.g. Hernquist vis-a-vis Sersic models. 
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Figure 6. Circular velocity profiles with radius, v c [r) = 
GM(r)/r, for cosmologically-motivated galaxy toy models. The 
stellar and halo masses are fixed to log 10 (M ic / Mq) = 10.85 and 
\og 10 (M v - lr / Mq) = 12.14, respectively, corresponding to egp = 
0.3. Three cases are considered with different adiabatic contrac- 
tion recipes, and are indicated by curves of different colours and 
linestyles, labelled by each recipe. For each case, a mass de- 
composition is shown, with the DM contribution as the bottom 
three curves, the stellar contribution as the intermediate steeply- 
declining blue dot-dashed curve (the same for all cases), and the 
total as the top three curves. A vertical dashed line shows the 
value of the effective radius adopted. 



including reasonable reproductions of the trend with M*.; 
there is no strong preference between the B+86 and G+04 
recipes. (The data clearly show a large variation around the 
models, part of which is just observational noise, and part of 
which is an important systematic variation that we will ex- 
amine in Section [4]) These conclusions are IMF dependent: 
with a Salpeter IMF, the /dm data points all shift down- 
wards to lower values, and the no-AC model is preferred 
on averagsQ (although there are suspicions of Salpeter not 
being a valid alternative: see Section [2. ip . 

For any fixed IMF, we see that it is a fairly generic 
expectation in a ACDM context for /dm to increase system- 
atically with Mi,, and to produce a DM-driven tilt in the 
FP. The reason is the aperture effect mentioned earlier: the 
scale-lengths of galaxies' stellar parts change more rapidly 
than for their haloes, causing the R e s region to encompass 

7 In paper I we focussed on the Salpeter IMF and stated erro- 
neously that the consistency of the data with our initial no-AC 
toy models was independent of IMF, when examining the relation 
between mass and mean DM density within 2 R s g. This relation 
for the data is indeed unaffected by IMF changes since both the 
implied masses and the densities vary with the IMF. However, 
we neglected to consider the effects on the models, where only 
the predicted densities change (because the R a g values at a fixed 
mass are changed: see further discussion in Section l3.2l l. 
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Figure 7. Central DM fraction versus effective radius. Black 
boxes with error bars show the average trend in bins, while smaller 
points show individual galaxies, colourised by the (smoothed) 
stellar mass according to the inset colour bar. /dm correlates 
with both M+ and R e g which are themselves correlated quanti- 
ties, but the tighter trend in this Figure compared with Fig. [2] 
indicates that i? cf f is the more important parameter driving the 
FP tilt. Sample ACDM toy models in bins of constant (as- 
suming egp = 0.1 and G+04 AC) are shown as solid curves: an 
increase of /dm with R c ff is naturally expec ted in these models. 
The d ashed curve shows the trend found bv lHumphrev fc Buotd 



increasingly large amounts of DM. This interpretation is 
supported by the fuvi-ReS trends in Fig. [Jj and by various 
FP analyse s in the literature based on s imula t ions of galaxy 
form at ion ( Bovlan-Kolchin et al] [2005. 200( 1; lOiiorbe et al.l 



20051, 12009; iRobertson et all l200rj ; Hopkins et al l |200S| : 



Covington uOod ). A si milar point from an observat ional ba- 



sis is made by Fig. 5 of iHumphrev fc Buotd (|201fl ). and we 
include their results for comparison in Fig. [JJ here the ab- 
solute values of /dm cannot be readily compared since their 
.Reff values come from the if-band and should be systemat- 
ically different from ours in the B-band, but their general 
trend with size seems comparable to ours. The role of _R c ff 
in more model-dependent DM scalings is considered further 
in Section \3. 21 

Systematic trends in esf can affect the slope of /dm- 
M* in either direction, but the effect is fairly weak. Even 
for a factor of 7 change in assumed halo mass, /dm changes 
only by ~ 0.1; i.e. large changes in overall halo mass do 
not translate to large changes in central DM content. To 
cancel out the DM-induced tilt, esf would have to system- 
atically increase with mass, which goes in the opposite di- 
rection to conventional findings, wherein the most massive 
systems have the highest virial M/L. If the FP tilt turns out 
to not be driven by /dm trends (as found in iTruiillo et al.l 
l2004l ; |Jun fc Imll200alAllanson et al.ll2009l) . then this would 



be inconsistent with vanilla ACDM expectations (i.e., where 
central DM densities scale slowly relative to the observed 
stellar densities). 



3.2 Dark matter profiles 

We next consider some interesting implications for DM den- 
sity profiles with radius. It is notoriously difficult to deter- 
mine the DM content in the centres of ETGs, and scant 
few studies to-date have been able to compare the empirical 
DM properties to ACDM theory. T+09 and paper I (Fig. 14) 
have found mean central densities for DM haloes that scale 
with mass roughly as expected for NFW haloes. Here we add 
the finding that the DM content also scales (to first order) 
with galaxy size as expected. 

We will see in Section 14.11 that much of the observed 
/DM-age anti-correlation can be traced to the i? c H-age anti- 
correlation — implying that our ACDM toy models are at 
least roughly correct. This point is illustrated more clearly 
in Fig. [8l where for one example mass bin, the stellar mass 
has been subtracted from the total mass, and the residual 
DM content is quantified as the mean value within R e s, 

By measuring the DM content of individual galaxies 
with similar masses at different radii, we are effectively able 
to map out a mean DM density profile over a factor of 3 in 
radius. This trick of studying composite profiles is formally 
valid only if the binned galaxies all have approximately the 
same DM profile — as might be the case for a no- AC NFW 
model. The condition is clearly violated for models that in- 
clude AC since the contraction varies among the differently- 
sized galaxies. However, the impact of this non-homology is 
reasonably small as discussed further below, such that we 
can get an approximate idea of the DM slopes, as well as 
more firmly test the null hypothesis of a universal NFW 
profilH 

The first point to notice from Fig. [8] is that the data 
for (pdm) and /dm (see Fig. [JJ have opposite trends versus 
R a B, which turns out to be crucial evidence for cuspy DM 
haloes on scales of ~ 2-15 kpc. As R e h increases for a fixed 
galaxy mass and halo profile, the amount of enclosed DM 
and therefore /dm both increase. On the other hand, the 
local DM density of a cuspy halo decreases with radius, so 
a larger R e g means a smaller (pdm) • 



8 Note this is different from the analysis of|Graham ct al. (2006b) 
who looked at (pdm) inside the effective radius of the DM rather 
than of the stars. Their derived slopes reflect the steeply declining 
outer regions of the DM haloes, while ours involve the central 
regions. 

9 One might suspect that at some level we are getting out what 
we are putting in, since our default dynamical models used in 
deriving /dm assume an a = 2 total density profile in order 
to extrapolate the (to measurements to i? c ff. However, we have 
confirmed that using the alternative constant-M/L profile yields 
similar results (a unchanged for Salpeter IMF, and steeper by 
Aa ~ 0.1-0.2 for Kroupa IMF). A galaxy sample with direct 
mass constraints nearer to R cfi (e.g. C+06; lAuger et alj [20091; 
lGravesll2009l) would be ideal for pursuing the density profile anal- 
ysis more robustly. We have also checked that these results are 
not sensitive to the stellar populations models (Appendix IB3L 
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Figure 8. Relation between DM density within 1 R c g and the 
-R c fj value for ETGs in one mass bin. The top panel shows 
the case of a Kroupa IMF, where the mass corresponds to 
log 10 (M* / Mq) ~ 11.25, and the bottom panel shows a Salpeter 
IMF for the same galaxies [log 10 (M*/M0) ~ 11.4]. Data for indi- 
vidual galaxies are shown as points coloured by their ages (after 
smoothing: see inset colour-bars). ACDM-based models are shown 
as curves: the lower blue curves are for no-AC, the middle black 
ones for G+04 AC, and the top red ones for B+86 AC. For each 
set of AC models, fixed values of esF = 0.7, 0.3, 0.1 are shown 
as solid, long-dashed and short-dashed curves respectively. Lines 
in the upper right of the top panel illustrate power-law slopes of 
a = — 1 and —2. The scatter in points (including the apparent 
plume to very low densities) could be entirely due to observational 
errors; see Fig. I All 



More quantitatively, we may consider a general power- 
law density scaling for the central DM halo, pdm(^) ~ r~ a . 
For a < 3, MdmO") ~ r s ~ a and therefore (pdm) ~ 
R cS ~ 3 M D M(r = Rcff) ~ RcH~ a - For an uncontracted NFW 
halo, a ~ 1 at the smallest radii, but near R c g in our toy 
models, a ~ 1.1—1.3 (stee per for higher masses because of 
their larger radii; see also lGraham et al"1l200 6a). AC effec- 
tively makes the DM even cuspier, with a ~ 1.6-1.8 for the 
G+04 recipe, and a ~ 1.8-1.9 for B+86 (see also Fig.©. The 
effect of constructing composite profiles from non-universal 
haloes is for the apparent slopes to be steeper by Aq ~ 0.3. 
This offset can be subtracted from the observational results 
to estimate the true slopes (see below), while more rigor- 
ous tests are made by constructing composite toy models 
for comparison to the data. 

Our observational results for a Kroupa IMF are illus- 
trated by Fig. [8] (top). The vast majority of the data points 
in this fixed mass bin do appear to clump around a com- 



mon DM mass profile, with a composite slope of a ~ 1.9 
(so we may infer a true slope of lvq ~ 1.6). A small minor- 
ity of galaxies may follow a shallower trend, with a residual 
trend of DM fraction versus age that is barely visible by 
the colourisation of the data points (to be discussed in Sec- 
tion The results for the other mass bins are similar, 
with a ~ 1.7-2.1, and suggestions of a shallower slope with 
increasing mass (see Fig. [9]), although it is hard to tell for 
sure since the exact slopes are sensitive to the details of the 
fitting procedur 

We also construct NFW-based model predictions as in 
Section ^. H and show them in Fig.JS] The steep slope for the 
majority of the points, along with the shallower slope with 
increasing mass, turns out to correspond nicely to a range 
of NFW+AC modelfEI with a minority of points consis- 
tent with uncontracted NFW haloes. The model interpreta- 
tions are somewhat degenerate to systematic variations of 
esF with R B ff (e.g. S+10 suggest an anti-correlation, which 
would imply the intrinsic slopes are steeper than they ap- 
pear). However, as seen in Fig. [8l the csf effects are fairly 
small, and overshadowed by our fitting uncertainties anyway. 

The {pom)-R c s relation is not a knock-down confirma- 
tion of ACDM but intriguingly it does seem to weigh against 
alternative models. If there is no DM, then of course we 
should have /dm = (pdm) = 0, although there might be a 
systematic error in the mass analysis, e.g. with the IMF. In 
this more general case we expect a ~ 3, and the apparent 
/dm to be constant with age and R e s — in strong contradic- 
tion to the data. 

The same argument may apply to, and pose a chal- 
lenge for, alternative gravity theories that seek to explain 
observational mass discrepanc ies at large galactic radii , 
usually in spiral galaxies (e.g., Brownstein fe Moffat! 120061 ; 
iFrigerio Martins fc Saluccil l2007h . It is beyond the scope 
of this paper to consider any of these theories in detail, 
but we will briefly comment on the most well s tudied case 
of M o dified Newtonian Dynamics (MOND; e.g. iTiret et all 
120071 ; ISanders fc Land! |200S| ; iFerreras et aT]|2009h . In this 
formalism, there is a characteristic acceleration scale ao ~ 
1.2 x 10 -8 cm s -2 , above which mass discrepancies should 
be weak or nonexistent. In the central regions of the ETGs 
studied here, the accelerations are large (~ 4ao) and the 
dynamical mass should agree well with the inferred stellar 
mass. 

Our observed mass discrepancies, along with the strik- 
ing R a f{ dependencies, would at first glance appear to rule 
out MOND without DM in ETGs, which might be re- 
lated to simil ar results in galaxy groups and cluster s (e.g . 



Sandersl2003l;|Pointecouteau fc Sirkl2005l ; lAngus et al.ll2008l ; 
Richtler et al.ll2008l ). However, we cannot make a definitive 



10 A fair amount of attention in the literature has been focussed 
on the slopes of total mass profiles rather t han the DM slopes dis- 
cusscd here (e.g. Ijiang fc Kochanekll2007l; iKoopmans et~aUl2009t 
iNipoti et a l. 2009; H umphrey 



BuotetoloT) . We have not made 



comparisons in detail, but our toy models imply total slopes at 
R e g of a ~ 1.8-2.2, becoming shallower with increasing mass. 
This trend is driven mostly by the stellar density becoming shal- 
lower, and is affected little by AC. 

11 The overall DM density shows a tendency to increase with 
mass, relative to fixed- esF models (see paper I), so the overall 
connection between AC and mass is not yet clear. 
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Figure 9. As Fig. [8] but with all ETG mass bins plotted to- 
gether (using a Kroupa IMF: see Fig[5]for colour key). In addition, 
data are shown for dwarf spheroid al galaxies (grey filled circles on 
the left; IWalker et al,|[20 09l. l201 0h and fo r late-type disk galaxies 
(black open circles: lMcGaugh et al.ll2007t V The ACDM toy model 
curves are as in Fig. [8] with the addition of sample grey and black 
model curves for the dwarf and disk galaxies. The dwarf models 
include a no-AC halo sequence with e SF = 0.04, 0.004, 0.0004, 
along with a case of esp = 0.004 with AC (G+04 recipe; solid 
curve). The spiral model uses a no-AC halo with egp = 0.3. 

statement without ca rrying out a careful MONDian analy- 
sis (see also Fig. 19 in lGerhard et alj|200ll . hereafter G+01; 
prelimin ary findings in | Romanowskvll200o1 : and forthcoming 
results in lCardone et al.ll2010l ). In particular, the systematic 
uncertainties in our Td yn and T* estimates must be taken 
into account, along with possible fine-tuning of the "inter- 
polation fu nction" between the Newtonian and MONDian 
regimes (cf. lMcGaugh|[20oi ). 

I f there are DM haloes w ith cored profiles (e.g. iBurkertJ 

ll995l ; ISaxton fc Ferr cras 2010), then /dm would still increase 
with Re,a, but (pdm) should be constant (a ~ 0). Again, 
using the wrong IMF could contribute an a ~ 3 effect, which 
in combination with a constant DM core could produce a 
shallower than a ~ 3 trend — but experimenting with IMF 
adjustments, we cannot recover a constant-core trend. 

Any of these scenarios could still be salvaged with epicy- 
cles, e.g. with systematic errors that correlate with galaxy 
size, but Ockham's razor suggests the first conclusion should 
be to prefer a standard cuspy DM model. Note that de- 
tailed dynamical models of individual galaxies have difficulty 
uniquely decomposing the stellar and DM mass profiles and 
thereby di stinguishing between co red and cusped DM haloes 
(cf. T+09; lNapolitano et ai]|2009l . hereafter N+09). Our ap- 
proach that includes stellar populations modelling in a large 
galaxy sample provides the first clear (albeit indirect) con- 
straint on the central DM profiles of ordinary ETGs. 

Before moving on, we consider the Salpeter IMF al- 



ternative, with data and models for the sample mass bin 
shown in Fig. [8] (bottom). The data are seen to typically 
lie in between the AC and no-AC models, with a slope of 
n ~ 1.7. Unlike the Kroupa case where both the /dm ampli- 
tude and the DM density slope match the AC models nicely, 
the Salpeter case shows less overall consistency with a simple 
ACDM model — although a caveat is that the frequent cases 
with apparent negative DM densities can make the results 
difficult to interpret. 

3.3 Context and comparisons 

There are several recent observational analyses of DM in 
galaxies that are relevant to our current study: we will make 
some comparisons first to studies of ETGs, and then to other 
galaxy types. We begin with two papers that are very sim- 
ilar in spirit to ours, as they combine estimates of total 
masses and stellar populations-based masses to derive cen- 
tral DM constraints in large samples of elliptical galaxies, 
and to compare these to ACDM-based toy models. 

T+10 studied a sample of strong gravitational lenses, 
using a combination of lensing and dynamics to derive the 
total masses. They found that a Salpeter IMF is generally 
consistent with a no-AC model, as did we. They did not 
consider AC models explicitly, and it seems very likely that 
(as we have found in this paper) lower-mass IMFs would be 
favoured if AC were included. 

These authors also found indications of systematic 
changes with galaxy mass that could be interpreted as a 
change either in IMF, or in the DM density slope such that 
it steepens with mass (a increasing). The latter effect as they 
mention could be due to AC becoming stronger with mass, 
and appears to go in the opposite direction of our results. It 
is beyond the scope of this paper to track down the reasons 
for this difference (their study did probe somewhat different 
regimes of radius and mass). 

S+10 derived their total central masses from stellar dy- 
namics, and used weak gravitational lensing to derive the 
total halo masses and concentrations (which we have had 
to assume in our model f^l) . Assuming a Kroupa IMF, they 
found relatively high central DM masses which they showed 
to be nicely consistent with G+04 AC models when breaking 
the sample down into several bins in mass and in size (the 
latter test being implicitly equivalent to our plots of (pdm) 
versus R c g as an alternative test of the DM profiles). Their 
large-radius constraints also allowed them to determine a 
correlation between size and virial mass at fixed K'h, which 
in the context of our formalism implies an anti-correlation 
between _R e ff and esF. If their sample has an i? e ff-age anti- 
correlation as in our sample, then this would imply an esf- 
age correlation. We will come back to this point in Sec- 
tion gT] 

In addition to the qualitative consistency of the S+10 
results with ours, their quantitative DM fractions of typ- 
ically /dm ~ 0.6-0.7 can be compared to ours over the 
same mass range [log 10 (M+/Mq) ~ 10.8-11.6] where we find 

12 The direct mass-concentration results of S+10 follow well the 
theor etical expectations for a WMAP5 cosmology l lMaccio et al.l 
2008), although one might expect ETGs to be systematically bi- 
ased toward higher-concentration haloes. 
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/dm ~ 0.4-0.6. This small difference might be accounted for 
by their sample selection of only round "central" galaxies — 
which are plausibly more DM-rich than flattened, satellite 
galaxies. These authors do note that their results seem to 
be somewhat high com pared to other literat ure studies with 
/dm ~ 0.3-0.5 (C+06: iGavazzi et al.l [20071 : see also G+01; 
T+09}. 

Ijiang fc Kochanekl (2007) analyzed a sample of galaxies 
using strong lensing and stellar dynamics, and ACDM mass 
models similar to ours. Using priors on either T* (equiv- 
alen t to a Kroupa IMF ) or total masses from weak lens- 
ing (jGavazzi et al.ll2007l ). they find that an AC model with 
esF ~ 0-3 is p r eferre d. The strong+weak lensing analysis of 
IGavazzi et all (|2007m on the other hand found consistency 
with a no- AC model. 

T+09 constructed detailed dynamical models of Coma 
cluster ETGs, where T* and Td yn were determined simulta- 
neously from the dynamics, assuming a parameterized DM 
model. They found relatively high DM densities (similar to 
ours), which they compared to cosmological semi-analytic 
models (using no- AC models for spirals as an additional con- 
straint, and adopting cored DM halo models for the ETGs 
which in this context tend to wnrfer-estimate the central DM 
densities). They concluded that B+86 AC haloes were pre- 
ferr ed over no- AC. 

iHumphrev et all |2009l) used X-ray emission from hot 
gas surrounding a small sample of massive ETGs to con- 
strain their mass profiles using ACDM-motivated models. 
With a Kroupa IMF, they found that no-AC models were 
preferred (although with somewhat high concentrations). 
One of their sample galaxies was also the subject of a de- 
tailed dynamical study using stars and globular clusters 
i|Shen fc Gebhardtj|2010h . whose mass results were system- 
atically higher and might very well be consistent with an 
AC model. This appears to be part of a significant unre- 
solved tension between dynamical and X-ray mass results. 
On the other hand, other X-ray studies fi nd surprisingly high 
DM concentratio ns for galaxy groups (|Buote et al.l 120071 : 
iDuffv et all 120081 ), which should generally overlap with the 
massive end of our ETC sample, and might support our 
findi ng that AC models are pr eferred. 

lLackner fc Ostrikerl ((2010) used toy models for "dissipa- 
tionless" and "dissipational" ETG formation which should 
be roughly equivalent to our no-AC and AC models. Adopt- 
ing the C+06 Y* values assuming a Kroupa IMF, they found 
that both types of models predicted too much DM in the 
galaxy centres. This disagrees with our conclusion that AC 
models predict about the correct amount of DM, and no- 
AC models not enough. Although there are some differences 
between our observational results and those of C+06 (see 
Appendix IB3|) . it may also be that Lackner & Ostriker did 
not allow for enough freedom in their halo masses and con- 
centrations to fit the data. Our own preliminary analysis 
of the C+06 data sug gested that a G+04 AC model would 
work well on average l|Roman owskv 200^). 

We next consider our ETG results in relation to other 
galaxy types, showing density-size relations for different 
ETG mass bins, along with data from other galaxy types 
in the literature, plotted together i n Fig. [9] We first of a ll 
consider 56 late- type galaxies from IMcGaugh et al.l l|2007l ). 
where dynamical masses were measured by classical rotation 
curves, and stellar disk masses by various methods includ- 



ing stellar populations modelling. We take results from the 
latter (with Kroupa IMF) and then infer the DM densities 
within R c s, which we define as 1.68 x i?d, where the disk 
sc ale-length Rd is t aken from various literature sources such 
as lMcGaughl l|2005l ). 

We plot these late-type galaxies in Fig. |5J and see that 
their DM haloes are ~ 5 times less dense than those of early- 
types at the same radii. Qualitatively similar conclusions 
were reached by G+01 and T+09, while we add the observa- 
tion that the DM density slopes are different: the late- types 
have shallower slopes than the early- types (a ~ — 1 versus 
a ~ —2). The ETGs could be brought into closer consis- 
tency with the late-types if a Salpeter IMF were adopted for 
the former, but if anything a reverse IMF trend might be 
expected a priori (as we will discuss further in Section f4. 41) . 

We also show in Fig. [5] a sample ACDM model curve 
for the late-types, for the case of log 10 (M* / Mq) — 10, 
£SF = 0.3 and no AC. This curve is very similar to the 
equivalent one for the lowest-mass bin of the ETGs, so for 
clarity we do not show additional model curves for the late- 
types (whose AC model curves should differ slightly from the 
ETG curves). The no-AC model appears to be preferred for 
the late-tvpej^|. echoing the conclusions of many other stud- 
ies dKassin et al.ll2006l; ibutton et~aj1l2007l ; IMcGaugh etTafl 
l2007l : lGnedin et alj|2007l ; IXue et aI1l2008l) — and contrasting 
with the inference of AC for the early-types when the same 
IMF is adopted. 

We next consider dwarf sphe roidals belong i ng to the 
Milky Way, from the compilat ion oflWalker et alj (|2009l ) (cf. 
the alternative compilation of I Wolf et al. 20ld ). This study 
extended upon earlier work th at found a common DM mass 
within a fixed physical radius (|Strigari et al.ll2~008l ) and sug- 
gested that these systems share a common "universal" DM 
profile with a — 1.6 ± 0.4 (see also iPenarrubia et al.|[201ol ). 
Plotting their data in Fig. we note first of all that the 
dwarf haloes do not appear to join up naturally with those 
of more massive galaxies (cf. McGaugh et al. 2010), and we 
infer that sampling a limited range of galaxy type, mass, 
and radius can readily produce the impression of a univer- 
sal profile. This conclusion is bolstered by our overplotting 
of a ACDM toy model for the dwarfs (using a fiducial stellar 
mass of log lo (M*/M0) = 5.5, Sersic index n — 0.25 and 
£SF = 0.004; a model with AC is very similar because of the 
DM dominance). For a factor of 10 variation in halo mass, 
the central DM densities change by only a factor of ~ 3, 
which is less than the scatter in the data. 

Finally, motivated by the work of iDonato et al.l (|2009l ) 
who claimed a universal D M projected density for all ga laxy 
types and masses (see also iKormendv fc Freemarj|2004l ). we 
calculate our own approximate DM surface density (£dm) 
from the data points in Fig. [9]by multiplying (pdm) by i? c ff , 
and showing the results versus stellar mass in Fig. QJj] Al- 
though our results are calculated within galaxy R e g rather 



13 However, IMcGaugh et al.l |2007) among others found that 
ACDM no-AC models do not match the data in detail, with 
cored DM halos being preferred. Note that if these haloes do have 
constant-density cores, then there is less concern about selecting 
the fairest radius fo r comparing their de nsities to ETG haloes. 

14 The (p) values in I Walker et all 11200911 appear to be total mass 
including baryonic contributions, but this should be only a ~ 10% 
effect. 
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discussion of the implications of our DM results for galaxy 
formation to Sections 14.31 and 14.41 
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Figure 10. Projected density of DM in the central regions 
of galaxies, versus stellar mass. The data and symbols are as 
in Fig. [9] For compa rison, the projected DM density result 
of iDona to et ahl 2003) is shown as a dashed line (see also 



I Gentile ^t**al*[l200E ; our density calculation is not exactly equiv- 
alent). The equivalent acceleration scale is provided on the right 
axis; the characteristic scale in MOND is log ao / (cm s — 2 ) ~ —7.9. 
Note that these quantities are all expressing the projection of the 
3D density within R e ff rather than a true total surface density, 
which would require model-dependent extrapolations outside the 
regions probed by the data. 



than the DM core radius as in Donato et al., we reproduce 
their results reasonably well, finding a fairly constant surface 
density on average across a large mass range of dwarfs and 
late types, but a lso suggesting a syste matic increase with 
mass as found bv lBovarskv et al] (|2009h . 

It is beyond the scope of this paper to consider ACDM 
theoretical predi c tions for the surfac e density in detail (cf . 

' 1 'Maccio et al l 120091: IKoposqv et al.1 [2009: 



Zha o et al] 120081: - _ _ 

Li et al.ll2009l;IOkamoto fc Frenkl2009l: Bovarskv et ahlbood : 



Cooper et al.ll2010l ; IStringer et al]|2010l ). However, the broad 



consistency of the uncontracted NFW models with the dwarf 
and spiral data in Fig. [9] suggests that the surface den- 
sity result is not particularly surprising within the standard 
paradigm, and does not necessarily imply cored DM haloes. 

We also find from Fig. \W\ that the ETG galaxies vio- 
late the constant density scenario for the other galaxies by 
a factor of ~ 10 on average, and a factor of ~ 5 in the same 
mass regime. This disagrees with the conclusions of Donato 
et al., and we note that their ETG results were based on a 
weak lensing analysis rather than central dynamical results 
equivalent to those used for the other galaxies. Other ETG 
dyna mical results from the literature support ours (G+01; 
T+09; iBovarskv et"al]|2009l; see also the strong lensing anal- 
ysis of 



Cardone. Angus, fc TortorgfolOl) . 



In summary, we can synthesize the constraints on the 
central DM content of all types of galaxies assuming a uni- 
versal Kroupa IMF, based on our study and on the litera- 
ture. We conclude that the early-type and late-type galaxies 
are broadly consistent with simple ACDM predictions with 
and without AC, respectively. 

The main exceptions to this emerging consensus are 
from X-ray studies (see above) and from the only com- 
pilation of large-radius dynamical results (N+09), which 
suggests very strong systematic AC differences among the 
ETGs. However, as discussed in Appendix IB2I selection ef- 
fects are suspected for the latter study. We postpone further 



4 INTERPRETING TRENDS WITH AGE 

At last we arrive at an analysis of the age trends. We provide 
an overview of the observations compared with basic ACDM 
toy models in Section 14.11 We then attempt to explain the 
age trends by the systematics of DM halo concentrations 
( Section [T2}, AC variations (Section 14. 3[) and non-universal 
IMFs (Section S3]). 



4.1 Overview of trends 

As previously mentioned, the variations in /dm at fixed 
appear to correlate with age, and comparison with the mod- 
els in Fig. [2] suggests that plausible systematic variations in 
esF with age would not be enough to account for the trends 
in the data. To see this more clearly, we construct models of 
/dm versus age in fixed mass bins, comparing these with the 
data in Fig. [TT] (keeping in mind that we expect some scat- 
ter in the data points from observational errors and from 
the necessary inclusion of a range of masses in each mass 
bin). 

Both data and models agree that /dm decreases with 
age at fixed mass, which can be seen as a natural conse- 
quence of the iieff-age trendj^l. However, the data in every 
mass bin show even steeper age trends than predicted by the 
models — particularly if we adopt models with AC, a process 
that couples the DM to the baryons and partially counter- 
acts the trend for small stellar sizes to produce small /dm 
(i.e. relative to uncontracted haloes, it leads to a shallower 
trend of /dm with age). The difference in slope between the 
data (d/oM/dage ~ —0.04 Gyr - ) and the simple models 
(~ -0.01 Gyr" 1 for AC, ~ -0.02 Gyr" 1 for no-AC) sug- 
gests a systematic variation with age of the DM itself. After 
correcting for the effects of R e g (in an inevitably model- 
dependent way), we support our initial conclusion from Sec- 
tion 12.21 that the central DM content is connected more 
closely to age than to mass, since /dm varies with age by 
up to ~ 0.3 relative to the model curves, and by only ~ 0.1 
with mass. 

We illustrate the effect of age in a different way in 
Fig. [12] by plotting the mean central DM density versus 
age, for both data and models in one mass bin. Because 
of the anti-correlation of R e g with age, the models predict 
a slight increase of (pdm) with age. However, the data show 
a decrease (also apparent from the colourisation in Fig. [8}, 
consistently with the results based on /dm- 

One obvious candidate to explain the DM-age trends 
would be esF, such that the younger galaxies simply have 
more massive DM haloes overall. This scenar io might be 
expected given cold-mode accretion at high-z (Deke l et al] 



5 IShankar fc Bernard I J2009l. hereafter SB09) also found a sim- 
ilar trend of excess dynamical mass with age, which they stated 
as driven mainly by the ii c ff variations. However, they did not 
demonstrate this conclusion explicitly, nor did they connect to 
specific DM models. 
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Age (Gyr) 

Figure 12. Mean dark matter density within 1 R e g as a function 
of stellar age, in the log 10 (M* / Mq) = 10.8 mass bin. Curves show 
model predictions as in Fig. 1111 while small points show empirical 
results for individual galaxies. Large points with error bars show 
the average trend in age bins. 
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Figure 11. Relation between central DM fraction and stel- 
lar age for ETGs in four mass bins: log lo (Af*/Af0) = 
11.6,11.2,10.8,10.3 from top to bottom. Data for four age bins 
are plotted as points with error bars as the mean and standard 
deviation for /dm- ACDM-based model predictions are shown as 
curves (see Fig. \8\ caption for details: the model trends with age 
are driven by R e f[ variations). The default Kroupa IMF is used 
here, but for comparison the results for a Salpeter IMF are in- 
cluded in the log 10 (M+ / Mq) = 11.2 panel: the large open boxes 
show the binned data values, and the lowest solid curve shows the 
egp = 0.7 no-AC model — which differs little from the equivalent 
Kroupa curve. Changes in the IMF affect the data much more 
than the models (where the IMF indirectly links the luminosity 



l2009al ). where star formation is thought to proceed very ef- 
ficiently; subsequent halo growth would have to include in- 
fall material with high esf in order to not overly dilute 
the net esp at low-z. However, although it is plausible that 
an esF-age correlation contributes to the /dm trendfj, we 
suspect it is probably not the dominant factor. From exam- 
ining Figs. [11] and 1121 systematic variations of a factor of 
10 or more in esp at fixed M* would be required. While 
there has been a great deal of attention given to associa- 
tions of galaxies with halo masses on average, there has been 
much less work on the scatter in the trends. It is outside the 
scope of this paper to investigate the plausibili ty of order- 
of-ma gnftude variations, but the new work of iMore et al.l 
|2010f) suggests a scatter in virial M/L of factors of only 
~ 2-5 for non-satellite ETGs (increasing with M*). 



4.2 Concentration variations 

We next consider possible connections between DM profiles 
and age, such that the amount of DM within i? c ff is not di- 
rectly coupled to the total halo mass. First of all, our NFW 
toy models assumed the exact same mass-concentration re- 
lation for all galaxy haloes. However, this should really be a 
mean trend whose intrinsic scatter correlates with the red- 
shift of halo collapse z c , since halo density basically reflects 
the background mass density of the universe at the time of 
collapse. 

This concentration-age systematic can immediately be 
seen to go in the wrong direction to explain the /DM-age 
data, if we make the reasonable assumption that stellar age 
correlates at least weakly with the assembly age of the DM 



16 S+10 found for a sample of low-z ETGs that halo mass cor- 
relates with R c ff at fixed stellar mass, by a factor of ~1.5— 2 in 
mass. If R e g anti-correlates with age in their sample as in ours, 
then egF correlates with age. 
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halo. Older galaxies should therefore have higher /dm and 
(Pdm) because of their denser, earlier-collapsing haloej^l. 

To illustrate this issue more quantitatively, we construct 
a modified version of our ACDM toy models. We begin with 
the simplest assumption that the age of a galaxy's stars is 
directly associated with its halo's z c . This is clearly not cor- 
rect in general since in the modern picture of "downsizing" , 
the star formation in massive galaxies largely precedes the 
DM asse mbly, and vice-versa in le ss massive systems (e.g. 
Fig. 5 of IConrov fc Wechslerj 1 2009 ) . However, it is beyond 
the scope of this paper to consider differential assembly his- 
tories at fixed final mass, so for now we will adopt our simple 
prescription as stated above. 

To a good approximation, haloes of all masses and at 
all times "collapse" with a fixed concentration of c = K ~ 4; 
the subsequent concentration evolution is basically mediated 
by the expansion of the universe. If a halo collapsed at z = 
2 C , then at z = its concentration will be c ~ K(l + z c ) 
with a scatter of Alog c ~ 0.05-0.06 (see IWechsler et al.l 
l2002l : lMacci6 et al.ll2008h . 

This line of argument implies that all galaxies with the 
same stellar age have the same halo concentrations, inde- 
pendent of mass. Certainly this is an oversimplified picture 
and might be undermined by late, minor gas-rich mergers 
biassing the inferred stellar ages. For now we will adopt this 
model as an ansatz and see where it leads us, with predic- 
tions illustrated by the model concentration-age bands in 

Fig.m 

To derive estimates from the data for concentrations, 
we proceed as follows. First, we collect the individual galaxy 
results into bins of age and stellar mass, calculating the av- 
erage /dm values (similar to Fig. II 1 1) . We then fit our ACDM 
toy models (Section [371]) to the /dm data. Given a fixed IMF 
and DM profile (NFW with or without AC), the remaining 
free parameters are now esf and c. We fix esf to be con- 
stant in each mass bin, and set it to a value that forces the 
average c (including all ages) to approximately agree with 
the theoretically predicted mean value. Note that for models 
with AC, the c value denotes the "original" concentration 
after correcting the observations for AC. 

The results of this fitting exercise for a Kroupa IMF 
and G+04 AC prescription are overplotted as data points 
in Fig. 1131 To reproduce the mean c-M v ir relation, we find 
£sf = 0.25 for the lowest mass bin, and esf = 0.1 for the 
others; it is noteworthy that such a simple, plausible ACDM- 
motivated model can be constructed to agree with both cen- 
tral and global DM constraints for ETGs. 

Now considering the age dependencies, we see that there 
is a clear trend for older galaxies to have less concentrated 
haloes — an effect that is diametrically opposed to our as- 
sumed trend for star formation to mimic halo collapse. If 
taken at face value, these concentration results would im- 
ply that ETGs on average formed their stars not long af- 
ter their haloes collapsed (z» ~ 1 versus z c ~ 2), but the 
late-collapsing haloes tended to form their stars very early 
(z* ~ 4 versus z c ~ 0.4) while the early-collapsing haloes 
formed their stars very late (z* ~ 0.4 versus z c ~ 5) — 



17 The stellar component could very well also be denser in older 
galaxies, but we are already accounting for this effect by including 
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Figure 13. The concentration of NFW haloes, as a function of 
halo mass. The underlying grey diagonal band shows the expected 
average trend for 68% of the haloes from ACDM simulations. The 
horizontal color bands show halo "ages" corresponding to collapse 
redshifts: 12, 11, 10, 8, 5 Gyr, or z c ~ 4, 3, 2, 1, 0.5, for red, 
pink, green, cyan, blue bands (top to bottom). The triangles with 
error bars show the data (assuming Kroupa IMF) in four mass 
bins, with age sub-bins of 12, 8, 4 Gyr (red, green, blue colours, 
respectively). Sample error bars are included for the intermediate- 
age bins, showing the scatter in masses and the 1 a uncertainties 
in the best-fit concentrations. These uncertainties mean that the 
one point with a very low concentration should not be taken too 
literally. 



although these latter trends are more apparent for the low- 
mass galaxies. We are not aware of any physical mechanism 
that would produce such counterintuitive variations at a 
fixed mass, and we suspect that concentration effects do not 
explain the observed DM-age trends. These conclusions are 
similar when adopting a Salpeter IMF with no AC. 



4.3 Adiabatic contraction 

The next alternative for systematic differences in the DM 
profiles is that the radial profiles have been systematically 
affected by baryons in different ways. The simplest picture 
in this context is for AC strength to vary with age. We do 
not explore this effect in detail, but refer to Fig. [TT] for a 
qualitative understanding. If the IMF has a slightly higher 
mass normalization than Kroupa (or if the esf values are 
somewhat lower than the range adopted so far), then the 
youngest galaxies are consistent with the B+86 AC recipe, 
with AC strength decreasing with age until the oldest galax- 
ies have had no AC. Variations in esf would then play a 
minor role in the central DM trends. In principle, we could 
confirm this interpretation by examining the (/9dm) radial 
profiles for different age bins, but the data are currently too 
noisy for this purpose. 

Is there any theoretical reason to expect AC to corre- 
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late with stellar age in this way? Or more generally, since 
we do not directly confirm the AC models, why should 
the DM be preferentially shifted away from the centres 
of "older" galaxies, and towards the centres of "younger" 
galaxies, by baryonic effects? One possibility concerns the 
smoothness of baryonic infall: the standard AC model is 
motivated by smooth gaseous dissipation, while clumpier 
collapse — both in the baryons and in the DM — is expected 
to diminish or even reverse the contra ction because of 
"feed b ack" from dynamical fric ti on (e.g. Debattista et al.l 
20081; iRomano-Dfaz et ail |200S| ; IJardel fc Sellwoodl 1200a 



Johansson et al.ll2009l ). 



Alternatively, more direct baryonic feedback could 
occur via outflows from AGNs and supernovae, poten- 
tially puffi n g up the DM halo (e.g . iBinney et al.l l200ll; 
Mo fc Mad I2OO4I; iMashchenko et all |200S|; jpeirani et ail 



20081 ; iGovernato et alj|20ld ; iDuffv et alj|201ol ). The notion 
of feedback effects on DM haloes is becoming commonplace, 
but the puzzle here is how to explain a systematic trend with 
time. Realistic simulations of galaxy formation including 
baryonic processes are only now becoming powerful enough 
to investigate such questions in detail. 

One specific scenario worth special mention is the 
emerging paradigm of high-z massive gala xy formation by 
cold streams of DM and baryons ( Dekel et al] l2009al lbl; 
lAgertz et aHl2009l ; iKeres et al.1l2009al lbh. The dumpiness of 
the infall in this picture could effectively bypass the contrac- 
tion process (as proposed bv lElmegreen et al.ll2008T ). These 
"wild disk" galaxies are though to initially form hot stellar 
disks and bulges in situ, and then evolve into low-z ellip- 
ticals, SOs, and early-type s pirals — whether by passive fad- 
ing o r by eventual merging (|Conrov et al.ll200Sl ; ICenel et all 
|200S| ). 

We would like to connect the DM-age trends for ETGs 
with the impacts of physical processes in the cosmologi- 
cal context of galaxy formation. To begin doing so, let us 
first consider even broader questions relating the different 
types of massive galaxies. As discussed in Section 13.31 our 
initial assessment of low-z galaxies is that the early-types 
have haloes with AC, and the late-types without AC. This 
type difference is not simply an age difference since among 
the ETGs, younger galaxies appear to have stronger AC. 
These observations raise an interesting cladistical question: 
if ETGs are generally formed in the mergers of spirals per 
the conventional wisdom, do they have roughly the same 
central DM densities? Also, if the most massive ETGs are 
assembled by dry mergers of smaller ETGs, are their respec- 
tive DM densities consistent with this picture? 

A related problem was raised long ago for the stellar- 
densities of galaxies. It is well known that in the merger of 
collisionless systems, the phase space density / cannot be 
increased, but observationally the central st ellar / values of 
ETGs are much highe r than those of spir als (|Carlberdll986l ; 
iHernquist et~aH 1 19931 : [Mao fe Mcl Il998h . This puzzle was 
solved by taking into account the dissipational infall of gas 
during a merge r, forming new stars in a high-density central 
starburst (e.g. [Kormendv fc Sanders! [l992l ; iBekki fc Shioval 
ll998l : ISpringej|200Ch . 

Returning to the DM issues, the 6D density / is of 
course not the same as the 3D density p, or 2D E, but we 



a 

S °8 



o 



1 1 1 r~ 



- 1 1 1 1— 



- 1 1 1 r~ 




•"Kit** 



.0 cfc o 



,0 



o 



■ o o 

0/ 

_l LiLl I I I I I I I ■* ■■' ■ I I I * I L 



10 11 
log 10 U./U B 



12 



Figure 14. Moan DM density within R c g, versus galaxy stellar 
mass. Filled dark circles are ETGs, while open bright circles are 
spirals, with colours indicating bins of constant size (i? c ft ~ 2, 
4, 9, 15 kpc, from top to bottom). Arrows indicate linear density 
growth with mass within the same radius. 



will not carry out a rigorous analysis of / herj"*l. Instead, 
noting that discrepancies involving / tend to be associated 
with p and E problems, we will take a very simplified ap- 
proach of conserving mass: the amount of central DM in a 
merger remnant should be rou ghly equal to the sum of the 
progenitors' central DM (e.g. iBovlan-Kolchin fc Mai |2004| ; 
iKazantzidis et ai1l2006l ). 

The observed DM density differences between early- 
and late-types appear from Figs. [5] and \TU\ to be factors of 
~ 5. However, in this context the figures are difficult to inter- 
pret because the two galaxy types occupy an almost disjoint 
region of Rcg-M* space. Comparison of the small samples 
of overlapping objects in this space indicate the real density 
differences may only be ~ 2 on average. 

We illustrate the trends in a different way in Fig. 1141 
where the mean densities are plotted against stellar masses, 
in bins of similar physical radius. The simplest expecta- 
tion here for galaxy mergers is that the DM density within 
the same radius grows linearly with mass, which the Figure 
shows is borne out remarkably well if considering the spirals 
as the progenitors of the ETGs. The full story is more com- 
plicated, as the true ETG progenitors are thought to have 
been high-z gas-rich systems with dense DM haloes (e.g. 
discussions in Section 14.21 and in G+01 and T+09). Also, 
there are potential effects in the mergers that could raise or 
lower the densities relative to the linear relation. A detailed 
theoretical exploration of DM density changes in mergers 



18 We have attempted rough estimates o f the coarse-grained 
phase -space density via {pdm)o"dm ( e -&- iDalcanton &; HoeanI 
1200 II) , but the results are very sensitive to how the DM dispersion 
i;dm is handled, and will require more careful analysis. 
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would be a useful contribution at this point, using also the 
constraint that the ETG haloes appear cuspier than those of 
the spirals. However, we can conclude that so far there does 
not appear to be an obvious problem in generating ETG 
halo densities from late-type progenitors. 

Now returning to the DM trends with age in the ETGs, 
progenitor bias seems an unlikely explanation. Higher DM 
densities in the early progenitors would produce the opposite 
low- z correlations to what we observe, and furthermore cur- 
rent observations of candidate high-z progenitors indicate 
relatively low centra l DM content, similar to low-z spirals 
ijBurkert et al.ll2010h . 

Instead, the possibility then remains that halo contrac- 
tion or expansion is a strong function of z. The scenario 
would be for halo contraction to become more important 
for ETGs that formed at later times, perhaps in a transi- 
tion from the cold stream phenomenon discussed above, to 
a more merger-dominated evolution that produces stronger 
contraction. Early, strong feedback might also play a role. 

An additional effect to consider is dissipationless merg- 
ing of ETGs: with repeated dry mergers, the initially seg- 
regated stellar and DM profiles should eventually converge 
to a well-mixed NFW-type profile for the total mass, im- 
plying a net migration of DM toward the galaxy's centre. 
One might surmise on this basis that older galaxies would 
have increased central DM content because of their longer 
post-star-formation per i od of merging. The simulations of 
iRuszkowski fe Springell ([2009) highlighted the effect of /dm 
augmentation through dry merging. However, the main /dm 
driver here may be the aperture effe ct due to the merger- 
induced Re tt growth (as the results of iBovlan-Kolchin et al.l 
120051 . 120061 suggest), rather than a significant evolution of 
the DM profile itself. 

4.4 Initial mass functions 

After exhausting a litany of possibilities for DM properties 
to vary with age, we consider the alternative that our DM 
inferences are wrong because we have assumed a universal 
IMF when estimating the stellar M/L values. The implica- 
tion would then be that the IMF varies systematically with 
stellar age, in the sense that younger galaxies have a higher- 
mass IMF: i.e., they have more stellar mass for a given lu- 
minosity, so the high dynamical M/L observations are due 
not to excess DM but to higher M*. In this context, when 
we couch the IMF in terms of "Salpeter" and "Kroupa" , we 
will not literally mean the same detailed IMF shape as those 
functions, but rather an equivalent overall normalization. 

The subject of IMF varia tions is highly contr oversial 
and unresolved: for a review see iBastian et ail (|201Ch . There 
are theoretical reasons to expect lower-mass I MFs for stellar 
populations that fo rmed in the early universe (|Larsonll2005l : 
iKlessen et al.|[2007l ). and observati onal suggestions for IMFs 
to be lower-mass at high er z (e.g. Ivan Dokkum|[200l : iDavel 
l200Sl ; lHolden et aUl201ch . These ideas would fit in well with 
the IMF-age trend we suggest above. All these results might 
further be consistent with suggesti ons that the IMF becomes 
more top-heavy at higher SFRs (IWeidner fc Kroupal 120061 : 
ICalura fc Mend 12009 : lHaas fc Anders! l2010h . since SFRs 
should have generally been higher at earlier times. Simi- 
larly, the more compact nature of the older ETGs might 
imply higher gas densities in the star-forming epoch, which 
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Figure 15. Stellar mass-to-light ratio produced by variable IMF 
(relative to Salpeter), versus stellar age. The small points show 
the results calculated for individual galaxies, the large points with 
error bars show median values in bins, and the dotted line shows 
the overall trend adopted. 



has also been suggested to produce more top-heavy IMF s 
jLee et al.ll2004 iMeurer et alJliooH : iKrumholz et alj|20ich . 

To be more quantitative about the implied IMF varia- 
tions, we construct a simple ad hoc model as follows. First 
we derive an average /dm-M* trend based on the data 
(Fig. [2]), and for each individual galaxy, solve for the Y* 
value that would agree with this trend (given its individual 
Tdyn value) . This value can be expressed relative to the ob- 
served Salpeter-based value by a fraction /imf = T*/T* >Sa i p 
(cf. T+10). The result of this exercise is shown in Fig. [15] we 
see a clear trend for the data to require a decreasing IMF 
mass with age, which we roughly parametrize by a linear 
function of /imp versus age. 

Given this new IMF assumption, we then re-derive some 
of our DM inferences. First we show the revised /DM-age 
trends in Fig. 1161 As intended by construction, the data 
show roughly constant DM fractions with age in fixed mass 
bins (and the galaxies with Salpeter IMFs are not in danger 
of being unphysical since these were the ones with appar- 
ently large /dm to begin with). As we have seen earlier, the 
/DM-age trends cannot be immediately interpreted without 
folding in the -Reg-age dependencies. Upon investigation, we 
find that with the varying IMF assumption, the strong -Reg- 
age anti-correlations (Fig. [5} have disappeared, so that there 
is on average no correlation. 

We find it a remarkable coincidence that the same IMF 
model removes the age trends in both /dm and -R G ff. This 
finding suggests that a systematic IMF variation could be 
real, and that it might account in part for the apparent 
tre nds with age and redshift of galaxy sizes (as suggested 
bv lMuzzin et al]|2009T ). Note that with our (somewhat arbi- 
trarily chosen) variable IMF, the /dm observations are now 
generally consistent with AC halo models. 

Given the opposing correlations of age with /imf and 
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Figure 16. As Fig. [4] with age-dependent IMF trend taken from 
Fig. OH 



M* (Fig. [3]), we might expect an anti-correlation between 
/imf and M*. We do see a weak suggestion of this ef- 
fect, which is a bit clearer when considering /imp versus <tq 
(Fig. I17|l . This result appears inconsistent with the strong 
positive /imf-co correlation found by T+10. However, these 
authors did not allow for AC in their models, which may 
impact their conclusions. 

As mentioned above, there are reasons to expect an 
IMF-age trend for ETGs qualitatively similar to our toy 
model interpretation. (Further analysis would be required to 
check for pro blems involving FP twisting with redshift; cf. 
Renzini 2006.) Now considering also late- types as relatively 
young, low-SFR, low-density systems, we would expect their 
IMFs to be more bottom-heavy (high /imf). Mergers of 
these systems to form the younger present-day ETGs would 
then also be expected to result in bottom-heavy IMFs, which 
is qualitatively consistent with the trend from Fig. 1151 but 
might be in conflict with IMF constraints in low-z galaxy 
disks. 



5 CONCLUSIONS 

We have continued an analysis begun in paper I of a large 
data-set of nearby early- type galaxies (ETGs), combining 
dynamics and stellar populations to constrain the central 
DM content. After having identified variations in DM as the 
main cause of the tilt of the fundamental plane, we have 
moved on to consider various scaling relations of the DM 
haloes, and connections to the star formation histories of 
the galaxies. 

Our basic observational findings are that the central 
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Figure 17. Stellar mass-to-light ratio produced by variable IMF 
(relative to Salpeter), versus central velocity dispersion. The large 
points with error bars show median values in bins, and the dotted 
curves shows the trend found by T+10, including their range of 
slopes. The overall normalization of our IMF is somewhat arbi- 
trary, so the key point of comparison is in the slope vs ctq. 



DM fraction /dm within an effective radius -R e ff has a strong 
anti-correlation with stellar age, and that the galaxy sizes 
also have an age anti-correlation. We have constructed com- 
posite profiles of DM density with radius, finding that they 
are on average cuspy, with inferred density exponents of 
~ —1.6 near R c ff. These profiles are steeper than literature 
findings for spiral galaxies, and the central DM densities of 
the early-types are denser overall, suggesting that gas-rich 
mergers would need to produce a net halo contraction. 

To further interpret the data, we have generated a series 
of ACDM toy models, including variable contributions from 
adiabatic contraction (AC). 

The results from comparisons of models to data are: 

• Models with AC fit well overall with a Kroupa IMF, 
while models without AC prefer a Salpeter IMF. 

• The size-age trends can explain part of the /DM-age 
trends: older galaxies show less evidence for DM because 
their more compact stellar centres probe less volume of the 
DM halo. 

• The remaining /DM-age trends are not easily explained 
by variations in halo mass or concentration, and suggest 
differences in baryonic effects on the DM, in the sense that 
younger galaxies have undergone AC while older galaxies 
have not. 

• An alternative scenario is for the IMF to be less massive 
for older stellar populations. 

There is ample scope for future insights and improve- 
ments. We plan to further investigate the galaxies' star 
formation histories in the context of theoretical mass as- 
sembly histories. Environmental trends can be investigated, 
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as these are expected to be important (e.g. Thomas et al 



2005 



2010 



Wechsler et ai1l2006l; iRogers et al.ll20ld ; iNierni et al 



La Barbera et al.l l2010bh \" Forthcoming high-quality. 



homogeneous, multiwavelength large surveys of low-redshift 
ETGs should also be able to refute, confirm or extend th e 
trends presented here (|Gravesll2009l : Icappellari et~alll2010h . 
Finally, the gold standard for probing galaxy mass profiles 
is extended kinematics data along with detailed dynamical 
modelling — both to provide more leverage on the DM inde- 
pendently of the stellar mass, and to sift indi vidual galaxies 
for the presence of a DM core or cusp (e.g. Thomas et al. 
20071: Ide Lorenzi et al.ll2009l : IWeiimans et al.ll2009l : iForestelll 



2009) 



Even if some of our current conclusions turn out 
to be completely wrong, we hope to have introduced a 
useful framework for interpreting mass results for large 
data sets of ETGs over cosmic time. The DM con- 
straints are part of a spectrum of clues that can ulti- 
mately be co mbined to pin down t he modes of ETG for- 
mation (e.g. ICovington et al.l 120081 ). Other avenues with 
consi derable promise include central rotati o n (Ijesseit et al.l 
extra light (Hopkins & HernauistJ l2010l ) and or- 



mi), 

bital st ructure (IBurkert et al. 20081 ) . as well as halo ro- 
tation dHoffman et alJ l2010l) and me tallicitv gradients 
ijWeiimans et al.ll2009l;lFoster et al.ll2009h. and globular clus- 
ter constra ints (|Rhode et all 120071 ; IShin fc Kawatal 120091 ; 
lBekkill2010h . 
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APPENDIX A: INVESTIGATING NEGATIVE 
DARK MATTER FRACTIONS 

Our derivation of /dm from estimates of Y* and Td yn yields 
a number of cases where galaxies have an unphysical /dm < 
0. These cases comprise 2%, 7% and 25% of the sample for 
the Chabrier, Kroupa and Salpeter IMFs, respectively. Here 
we investigate whether these fractions could be compatible 
with simple observational scatter in and Td yn - 

We have estimated the uncertainties in to be ~ 15% 
using Monte Carlo simulations of SED fitting. For Td yn , we 
may generically consider it as derived via a virial relation: 
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where K is a virial coefficient, G is the gravitational con- 
stant, and L is the luminosity (see C+06 equation 19). In 
this context, the value of R e g is arbitrary, and the uncer- 
tainty in the physical units R E g/L is dominated by the 
~ 10% distance uncertainty. To estimate cr e ff we first con- 
sider the central <to with a measurement error of ~ 5%, and 
extrapolate this to the value averaged over R c g — an exercise 
with a ~ 5% uncertainty (C+06; we consider systematics 
separately). The value of K was found by C+06 to have a 
~ 15% scatter. The net uncertainty in Td yn is then ~ 23%, 
and the uncertainty on /dm varies from ~ 0.1 to ~ 0.2 for 
high to low /dm values. 

We begin by assuming some fraction of the galaxy sam- 
ple have true /dm ~ 0.1, which we find is the lowest plausible 
value in a ACDM context fSection l3.1|) . Using the uncertain- 
ties discussed above, we generate a random sample of /dm 
measurements and scale a histogram of their frequency to 
reproduce the number of /dm < observations while not 
exceeding any of the /dm > frequencies (Fig. IA1|) . For 
the Chabrier and Kroupa IMFs, this scaling process implies 
that ~ 15% of the sample have /dm ~ 0.1. For the Salpeter 
IMF on the other hand, it is difficult to find a scaling that 
reproduces the negative- /dm tail without violating the ob- 
served distribution to higher values; the best fit has ~ 55% 
with /dm ~ 0.1. 

This difficulty with Salpeter would be eased if a true 
/dm = is assumed for half the galaxies, if there were sys- 
tematic errors at the level of ~ 15%, or if the random T* 
measurement errors were actually at the level of ~ 25%, or 
were non-Gaussian. We would not at present rule out any 
of these possibilities and so cannot categorically exclude a 




Figure Al. Distributions of DM fraction for early- type galax- 
ies. The open histograms are the observational results, with a 
different IMF for each panel (as labelled). The shaded green his- 
tograms are random realizations of galaxy subsamplcs having in- 
trinsic /dm = 0.1, with error bars illustrating the approximate 
Poissonian uncertainties. See text for further details. 



universal Salpeter IMF. Recalibrating our Td yn and T* esti- 
mates to external results (as discussed in paper I Appendix 
A) would on the other hand make the situation more prob- 
lematic for Salpeter by increasing the number of objects with 
/dm < 0. 



APPENDIX B: CROSS-CHECKS ON AGE 
DEPENDENCIES 

Here we carry out various tests on the robustness of the DM- 
age trends found in Section T2.2I In Section IBll we compare 
our /dm results to other literature results, and in Section lB2l 
we examine the implications of results at larger radii. We 
explore systematic effects in our stellar populations models 
in Section TB3I 



Bl Central DM content 

We check here whether our results on central /dm and age 
are consistent with other results in the literature. First we 
consider the SDSS-based analysis of SB09. They estimated 
central dynamical and stellar masses for a large sample of 
ETGs at z ~ 0.1-0.2 using somewhat different techniques 
to ours in paper I. After correcting the galaxy luminosities 
to a common age, they found total M/L trends that depend 
sy stematica l ly on age (their Fig. 2). Using the information 
m iBernardil (2009) to convert these corrected luminosities to 
stellar masses, we convert their M/L values to Md yn /Af* and 
then to /dm- It also appears that they assumed a Chabrier 
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Figure Bl. DM fraction with age of z ~ 0.1-0.2 ETGs, in bins 
of stellar mass, from SB09; the data have been converted from 
using a Chabrier IMF to a Kroupa IMF. The colours show the 
same mass bins as in Fig. 1111 For comparison, the lines show our 
median results in each mass bin. Also shown as large stars are 
typical results from low - and high-z ETG samples from C+06 and 
ICappellari et al.l (|2009h . The "ages" plotted in this case actually 
correspond to the redshifts of the observations, but we do not 
have enough information to correct these to z = ages (which 
will qualitatively be higher than the plotted points). 



IMF, so we further convert their results to a Kroupa IMF 
for comparison with our default models and results. 

The results of this exercise are plotted in Fig. IB1I there 
are indications of an anti-correlation between /dm and age, 
particularly for the higher mass bins. For comparison, our 
results are also plotted; other than a 30% overall offset in 
M/L whose origin is unclear, the results are generally con- 
sistent. The SB09 data do suggest somewhat shallower cor- 
relations that might be explained by a pure R e g-age effect 
(as stated by these authors) but they do not cover a large 
enoug h range in age to b e sure. 

iGraves et alj (|2009bl ) carried out a different analysis of 
quiescent ETGs in the SDSS. They again analyzed the cen- 
tral dynamical and stellar masses separately, and mapped 
various stellar populations parameters on approximate slices 
of the FP defined by ao and R c g. Although these results do 
imply a systematic correlation between size and age at fixed 
mass, if one considers a single ao-Raff grid-point, then one 
can control for size and mass dependencies and look for any 
residual correlations in the perpendicular direction of sur- 
face brightness. 

These authors found a negative correlation between age 
and surface brightness, which because dynamical mass is ap- 
proximately constant implies that there is a positive correla- 
tion between age and overall M/L. They mentioned that in 
a forthcoming paper they will find that stellar populations 
effects are not enough to explain this age-M/L correlation, 
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Figure B2. Our results for central DM fraction versus age, 
bins of constant <r (~ 120, 180, 250, 300 km s" 1 ). 



which therefore implies a net positive age- /dm correlation. 
This result is opposite to the trend that we find for a residual 
negative age- /dm correlation after the age-i? c ff correlations 
are accounted for. It remains to be seen exactly how and 
why our results differ. 

iBastian et all (|20T0t ) Fig. 4 provides an intriguing com- 
pilation of comparisons between T d vn and T* at z = an d 
z ~ 2 for ETGs from C+06 and ICappellari et al.l |2009h . 
While there are many issues in making such a direct com- 
parison between the two galaxy samples the apparent trend 
does agree qualitatively with our results, although we can- 
not evaluate the potential residuals from a size-driven trend 

(Fig. [mi 

since Graves et all (|2009bl ) among others emphasize the 
primary importance of ao (rather than K'h) in correlating 
with stellar populations parameters, we show in Fig. IB2l our 
results for /oM-age in bins of constant ao. The trends are 
similar to those found in bins of constant M*. However, we 
do not use ao in general in this paper because it is less 
straightforward to incorporate in the toy models, and be- 
cause /dm depends on a dynamical mass determination that 
is intrinsically correlated with ao. 



B2 Large radius DM content 

Next we consider DM conclusions from dynamical studies 
extending to large galactocentric radii, quantified as the 
M/L gradient parameter introduced in iNapolitano et all 
(2005). The data-set of 25 galaxies is from paper I (Fig. C2, 
with an update on NGC 4374 from Napolitano et al., in 
prep) , where we confirmed that the gradient correlated well 
with /dm from the central regions. We show some of our 
toy model predictions from this paper in Fig. IB3I (upper left 
panel). In the absence of any residual DM-age trends, the 
gradient is expected to decrease with age, since it is defined 
relative to the galaxy R c s, which in turn decreases with age. 

At first glance, the data appear to support this ex- 
pectation, with the residuals unclear because of the scat- 
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ter. However, closer examination of the data-sets plotted 
raises a red flag: the four "young", DM-dominated galax- 
ies are all well-known nearby massive group- and cluster- 
central ellipticals that are normally thought to have very 
old stellar populations (averaged over R c ff), not with ages 
~ 4-5 Gyr (NGC 1407 M49 M87, NGC 5846; e .g. C+06; 
ISanchez-Blazquez et aill2007l ; ISpolaor et afll200Sf ) . 

It turns out that these galaxies were all fitted with an 
unrealistic super-solar metallicity [Z = 0.05) which is prob- 
ably a reflection of the fundamental age- metallicity degen- 
eracy in stellar populations analyses. We re-ran our anal- 
ysis with all the metallicities fixed to Z — 0.02, which for 
the four problematic galaxies yielded more credible ages of 
14 Gvil 19 !. The revised M/L-gradient results are shown in 
the top right panel of Fig. IB3I where it now appears that 
the DM content may increase with age. The first thing to 
keep in mind when considering this apparent inconsistency 
with our main /DM-age result is that the DM content within 
-Reft and within ~ 5 -R ff may very well not be tightly corre- 
lated. The central DM content may be less a reflection of the 
overall DM content and more closely related to the details 
of the baryonic-DM interplay at the centres of haloes. 

The main effect however appears to be small number 
statistics coupled with selection effects. In the bottom pan- 
els of Fig. IB3I we show the /dm results for both metallicity 
assumptions. The stellar M/L is not very sensitive to the 
age-metallicity degeneracy, and in fact is affected in a way 
that roughly parallels the overall /DM-age trend: higher age 
and lower Z yield higher M/L* and lower /dm- Thus the 
general trend for our full sample is not qualitatively affected 
by changing Z, but quantitatively shifts to larger ages over- 
all (see also Fig. IB4|). 

The subsample of galaxies that have large-radius data 
is not so fortunate. When fixing Z = 0.02, their /DM-age 
results buck the overall trend, just as seen from the large- 
radius results. This sample has only one "young" object (<^ 7 
Gyr), and the old objects are dominated by systems like 
M87 which are known to reside at the centres of massive 
groups and clusters, and appear to comprise the high- /dm 
tail of the overall distribution. Such galaxies are some of 
the first targets of large-radius dynamical studies because 
of their rich supply of mass tracers such as globular clusters 
and planetary nebulae. The price to pay from selecting such 
systems is that they could provide a very biassed view of the 
Universe — as now appears to be the case when considering 
central DM content. 

This selection effect may explain the curious DM halo 
mass-concentration trends found by N+09, and demon- 
strates the importance of constructing an unbiassed galaxy 
sample. For now, the DM-age implications from large-radius 
tracers are totally inconclusive, and will require completion 
of a large systematic survey to make any progress. In the 
meantime, valuable spot-checks could be provided by study- 
ing the large-radius DM content of a few galaxies that help 



drive the apparent /oM-age relation, e.g. young DM-rich sys- 
tems like NGC 3626 or old DM-poor ones like NGC 7454. 



9 For galaxies with a previous best fit of Z = 0.02, the age and 
T* results with fixed Z are not exactly the same as before. This 
is because our procedure does not involve a simple best fit, but 
rather a Monte Carlo approach using the median of a distribution 
of best-fit values: see paper I. 




Figure B3. DM-age trends in a subsample of galaxies with large-radius dynamical tracers. The top row shows the "mass-to-light ratio 
radial gradient", which as an observational parameter is based on dynamics only and is independent of the IMF. Colourised model curves 
and mass bins are as in Fig. 1111 The models are for the case of Kroupa IMF using a G+04 AC recipe, with line-styles showing different 
£SF values (note that AC has only a small impact on these model predictions since it is a process that acts strongly at small radii). The 
bottom row shows the DM fraction within R e ff , with the small symbols for our full galaxy sample (compare Fig. |4j and large symbols 
for the subsample with large-radius results. The left panels show results based on our general stellar populations models with metallicity 
left as a free parameter, and the right panels show the results with metallicity fixed to Solar (Z = 0.02). 
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Figure B4. Trend of DM fraction with age, for different stel- 
lar populations assumptions. The solid black line is our stan- 
dard model, and see Fig. A2 of paper I for explanations of the 
other line sty l es. An additional model not included here is from 
IConrov et al. 1 1120091 1, but for the wavelengths and ages used here, 
this would be equivalent to BC03 (see their Fig. fO). We have also 
tested delayed-exponcntial and truncated (step-function) SFH 
models, which are not shown here but again do not substantially 
change the trends above. 



B3 Modelling systematics 

We next consider whether systematic uncertainties in our 
stellar populations modelling could be affecting our DM in- 
ferences. We have already explored this issue for the /dm- 
mass relations in paper I; here we consider /oM-age. Age 
and stellar M/L are both derived from the same models, 
and a positive error in age would correlate with a positive 
error in stellar M/L which would produce a negative error 
in /dm and thus mimic an anti-correlated /oM-age trend 
solely because of the correlated errors. 

To see if such a systematic could account for our ob- 
servational result, we experiment with different stellar pop- 
ulations assumptions, showing the results in Fig. IB4I The 
only cases where the /oM-age trend is appreciably different 
from our standard estimate is when we unrealistically fix 
Z = 0.05 (green lines). 

Given the problem identified in Section IB2I where al- 
lowing the metallicity to vary freely can result in too many 
"young" galaxies, we consider in particular our results when 
we fix Z — Zg. The implications for /dm were shown in 
Figs. IB3I (right panels) and IB4I (solid orange line), and we 
also show the overall (pDM)-fleff trend in Fig. IB5I The re- 
sults are similar to our standard model. 

As a reminder, in paper I (Appendix A) we compared 
independent estimates of Td yn and to C+06 for galaxies 
in common. Our T* values for the same IMF were ~ 20% 
lower, and our Td yn values were higher (by ~ 10% for the 
brighter galaxies and ~ 30% for the fainter galaxies). Thus 
our /dm values are systematically higher, becoming more 
discrepant for fainter galaxies (up to A/dm ~ 0.25). The 
C+06 results on their own may imply consistency with no- 
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Figure B5. Central DM density versus effective radius for our full 
ETG galaxy sample. Red points show results using our standard 
assumptions in the stellar populations models, and blue shows 
the case where metallicity is fixed to solar. 



AC models fo r the fainter galaxi es, and strong AC for the 
brighter ones (|Romanowskvll2006l ) . which would be in better 
agreement with the results of N+09. We have also checked 
the /DM-age trends using a recalibration as discussed in 
paper I, and found that for the fainter galaxies, the anti- 
correlation becomes somewhat steeper, so our basic /DM-age 
result does not go away. 

An issue not discussed in paper I is the potential effect 
of AGN emission on the observed galaxy colours and thus on 
the T* and age inferences. The AGN colours would gener- 
ally mimic a stellar population of ~ 3 Gyr age and skew the 
inferences toward younger ages and higher /dm- Matching 
the SPSS DR4 ETG ca talogue with the AGN catalogue of 
iKauffmann et all (|2003f ). we find that ~ 15% of the galax- 
ies in our sample's mass range have s trong AGNs (Type I I 
with strong [O III] emission; see also iGraves et ai]|2009a|) . 
However, based on the AGN study of lSchmitt et all l| 19991) . 
we estimate that the contaminant light would account for 
only ~ 1-2% of the total light within R e g , and would bias 
our results by at most 5%. 



